Non-ergodic Intensity Correlation Functions for Blinking Nano Crystals 
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We investigate the non-ergodic properties of blinking nano-crystals using a stochastic approach. 
We calculate the distribution functions of the time averaged intensity correlation function and show 
that these distributions are not delta peaked on the ensemble average correlation function values; 
instead they are W or U shaped. Beyond blinking nano-crystals our results describe non-ergodicity in 
systems stochastically modeled using the Levy walk framework for anomalous diffusion, for example 
certain types of chaotic dynamics, currents in ion-channel, and single spin dynamics to name a few. 
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Statistics of fluorescence intensity signals from single 
molecules, atoms, and nanocrystals are in many cases an- 
alyzed using intensity correlation functions (e.g., QQ])- 
These correlation functions are used to investigate a 
wide range of dynamical behaviors, for example effects 
like anti-bunching and the stochastic dynamics of large 
molecules coupled to their environment (0, |j| and Ref. 
therein) . In standard theories it is assumed that the pro- 
cess of photon emission is stationary and ergodic. In 
contrast, measurements of intensity correlation functions 
obtained from single nano-crystals (NCs) exhibit a non- 
stationary and non-ergodic behavior y, |y| ; as such these 
systems exhibit a statistical behavior very different than 
other single emitting objects. 

More speciflcally, the fluorescence emission of single 
colloidal NCs, e.g. CdSe quantum dots, exhibits inter- 
esting intermittency behavior Q- A standard method to 
analyze such intensity signals is to define a threshold Ith 
and define tvi^o states: on if I{t) > Ith and off otherwise. 
For capped NCs (e.g. CdSe(ZnS) core-shell NC) on and 
off times exhibit power law statistics [slal; their proba- 
bility density function (PDF) behave like tpiTLL ot r~(^+^^ 
for large r, and 6 < 1. For example in ja| 215 NCs 
were measured and the exponents 9on = 0.58 ± 0.17 and 
9off = 0.48 ± 0.15 were found (note that within error 
of measurement 9on — Ooff — 0), further all NCs are 
reported to be statistically identical. Since 9 < 1 the 
average on and off times are infinite. The divergence of 
occupation times naturally leads to non-ergodicity in the 
blinking NCs ^, ^j (see also [l^l for related discussion) . 
Other measurements classify the intermittency based on 
time average correlation function or closely related power 
spectrum [5, 11, 12|. In this type of analysis of the ex- 
perimental data there is no need to introduce a threshold 
value Ith- 

From a single realization of intensity trajectory I{t), 
recorded in a time interval (0,T'), we may construct the 
time averaged (TA) correlation function 



CTA{t',T') 



lolit^ 



t')I{t)At 



T 



(1) 



where we denoted T = T' — t' . On the other hand we 
may generate many intensity trajectories one at a time, 
and then average to obtain the ensemble average corre- 
lation function {I(t)I{t + t')). Single molecule experi- 



ments investigate time average correlation functions. If 
the random process I{t) is ergodic, the time average and 
the ensemble average correlation functions are identical 
in statistical sense, provided that the measurement time 
is long. Theories of correlation functions of single emit- 
ting objects are many times based on the assumption 
that the single molecule intensity trajectories are ergodic 
and hence for simplicity, theories concentrate on the cal- 
culation of the ensemble average correlation function. In 
this Letter we quantify the non-ergodic properties of the 
correlation functions of bhnking NCs, using a stochas- 
tic approach. A related question of non-stastionarity, or 
more specifically aging of the ensemble average correla- 
tion function is a subject of intensive research in the 
literature Hill 111 111. 



We use a simple two state stochastic model, with which 
correlation functions and non-ergodicity of the NCs are 
investigated. The intensity I{t) jumps between two states 
I{t) — 1 and I{t) — 0. At start of the measurement t = 
the NC begins in state on /(O) = 1. A schematic re- 
alization of the intensity fiuctuations is shown in Fig. 
ni The process is characterized based on the sequence 



{rf 



off 



off 



■} of on and off sojourn times 



or equivalently according to the dots on the time axis 
ti,t2,---, on which transitions from on to off or vice 
versa occur (See Fig. QJ. The sojourn time t^ is an 
off time if i is even, it is an on time if i is odd. The 
times Ti are drawn at random from the PDF iP{t). These 
sojourn times are mutually independent, identically dis- 
tributed random variables. We will consider the case 
V'(r) ex T-(i+«) where < 6* < 1. When 6* > 1, or 
when iP{t) is exponential, we find an ergodic behav- 
ior. Note that similar intermittency behavior describes a 
wide range of physical systems and models. In particular 
the model describes a Levy walk, which is an important 
stochastic model for anomalous diffusion describing dy- 
namics of systems mentioned in the abstract (2. Iiall7l|. 
The power law behavior of the sojourn times of the NCs 
can be explained based on simple physical models. For 
example random energy trapping models and random 
walks models [l^ can be used to model such behavior, 
though currently it is not clear yet what is the physi- 
cal mechanism behind the observed power law behavior 
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In Fig. m ten typical simulated correlation functions 
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Figure 1: (a) A schematic representation of intensity blinking, 
note that we redefine in to be equal to T. (b) Ten numerically 
generated realizations of the correlation function CTA{t' ,T') 
versus r — t' /T' , for 6 = 0.8. The correlation functions ex- 
hibit non-ergodic behavior and are random, for ergodic pro- 
cesses all the ten time averaged correlation functions would 
follow the same master curve, namely the ensemble average 
correlation function. 



are plotted, the most striking feature of the figure is 
that the correlation functions are random. These cor- 
relation functions are similar to those obtained in the 
experiment Q. Mathematically, the question of non- 
ergodicity may be formulated in the following way. Since 
the process I{t) is random the time average correlation 
function CTA{t',T') is random. For ergodic processes, 
and in the long measurement time limit, the distribution 
of CTA{t\T') is delta peaked and centered around the 
ensemble average correlation function. For non-ergodic 
processes the goal is to obtain the non-trivial limiting dis- 
tributions of CxA^t' ,T') which differ from the narrowly 
peaked delta functions found for ergodic processes. In 
what follows we will denote PcTA{t' T')ix) to be the PDF 

oiCTA{t',r). 

To start our analysis we rewrite the time average cor- 
relation function as 



CTA{t',T') 



T 



(2) 



where we used the initial condition that I{t) = 1 at time 
t — 0. Hence I{t) = 1 in ii_i < t < ti when i is odd, 
otherwise it is zero. The summation in Eq. Q is over 
odd z's, and i„ = T, namely n— 1 in Eq. IpJ is the random 
number of transitions in the interval [0,T]. From Eq. (01 
we see that the time averaged correlation function is a 
sum of the random variables 



n-t' + lit' i odd U - ti-i > t' 
I{t + t')dt = { liTi i odd U - U-i < t' 

i even 



/;-+*' iwdt 



i,.^ 



f''+'' i{t)dt 



if U - ti^i > t' 



iiU-ti^i < t'. 



(4) 



The li's are time averages of the signal I{t) over periods 
of length t' or n —ti — ti-i. Using Eqs. I|2I3|I we find an 
exact expression for the correlation functions in terms of 
{n} and {li} 



TCTA(t',T') 



>>- 


^ ^ (l-J.)r,-i' Y. 


odd 


i odd i odd 




n <t' n> t' 



a -I,' 



(5) 



(3) 



The first term on the right hand side of this equation is 
T+ the total time spent in state on in the time interval 
[0,T], in the remaining two terms we have considered 
sojourn times Ti larger or smaller than t' separately. 

We now illustrate the rich behaviors of the PDF 
PcTAit',T'){x) using numerical simulations, and later we 
consider the problem analytically. We generate random 
realization of the process using tp{T) = 9t~^~^ for t > 1 
and show two cases: — 0.3 in Fig. Hand 9 — 0.8 in Fig. 
El In both figures we vary r = t' /T' . The diamonds are 
numerical results which agree very well with the theoret- 
ical curves, without any fitting. First consider the case 
r = in Figs. I2I3I For 9 = 0.3 and r = we see from 
Fig. El that the PDF PcTA(t',T'){z) has a U shape. This 
is a strong non-ergodic behavior, since the PDF does not 
peak on the ensemble averaged value of the correlation 
function which is 1/2 for this case. On the other hand, 
when 9 = 0.8 the PDF Pc^Ait' ,t'){z) has a W shape, a 
weak non ergodic behavior. To understand the origin of 
this type of transition note that as ^ we expect the 
process to be in an on state or an off state for the whole 
duration of the measurement, hence in that case the PDF 
of the correlation function will peak on CxAit' ,T') — 1 
and CTA{t' , T') — (i.e U shape behavior). On the other 
hand when -^ 1 we expect a more ergodic behavior, 
since for 9 = \ the mean on and off periods are finite, 
this manifests itself in a peak of the distribution function 
of CTA{t', T') on the ensemble average value of 1/2 and a 
W shape PDF emerges. Note that for 9 < 1 there is still 
statistical weight for trajectories which are on or off for 
periods which are of the order of the measurement time 
T', hence the distribution of Cta(0, T') attains its max- 
imum on Cta{0,T') = 1 and Cta(0, T') = 0. For r > 
we observe in both figures a non-symmetrical shape of the 
PDF of the correlation function, which will be explained 
later. 

We first consider the non-ergodic properties of the cor- 
relation function for the case t' — 0. It is useful to define 

I[a,b]= I imt/ib-a), (6) 

J a 

the time average intensity between time a and time b > a. 
Using Eq. Q and for t' = the time averaged correlation 




Figure 2: The PDF of CTA{t' ,T') for 9 = 0.3 and different 
values of r = t' /T' . The diamonds are numerical simulations 
and the curves are analytical expression obtained for: (a) 
r = 0, Eq. © solid curve, (b) r = 0.01,0.1, Eq. lO dashed 
curve, and (c) for r = 0.5, 0.9, 0.99, Eq. I|15|l solid curve. In 
the ergodic phase the PDF of CTA{t' ,T') would be peaked 
around the ensemble average correlation function, which for 
r = falls on 1/2 and for i' ^ cxj is on 1/4 (for any r / 0). 
We see that any measurement is highly unlikely to yield the 
ensemble average when = 0.3. 



function is identical to the time average intensity 



Cta{0,T) 



'-[0,T] 



11 



(7) 



The random correlation function Cta(0, T) has a known 
asymptotic behavior in the limit T ^ oo, found originally 
by Lamperti [l^ (see also [IjI); this PDF is denoted with 

limT^oo PcTAio.T){x) = hix), and 



We now consider an analytical approach for the case 
t' < T. The behavior of PcTA{t',T)l^) for i' ^ is non- 
trivial, because the li's in Eq. ^ depend statistically 
on the random variables r^ . To treat the problem we use 
a non-ergodic mean field approximation. We noticed al- 
ready that Ti defined in Eq. Q are short time averages 
of the intensity, hence using mean field theory approach 
we replace the Ji in Eq. Q with the time average inten- 
sity 1(0, T] ; specific for a given realization. Replacing Xi 
with the ensemble average intensity is not appropriate. 
Hence within mean field 



TCTA{t',T')^I[0,T]T^{l-I[0,T]) {t'N+ 



E 

i odd 
n < t' 



(9) 
where iV+ is number of odd (i.e. on) intervals satisfying 
Ti > t' and i < n. 

We now investigate the distribution of CxAit' ,T') us- 
ing the approximation Eq. Q, leaving certain details of 
our derivation to a longer publication. First we replace 
TV"*" with its scaling form. Let P{t > t') — J^, V(''')dT 
be the probability of r being larger than t' , we have 

N+ ~ KP{t > t')T+/ J^'^ TiP{T)d T, where if is a con- 

stant of order 1, and T^ / /g Tip{T)d t is total number 
of jumps in time interval T'^ . A more refined treatment 
yields 



N^ 



sinTT^ 
Tr9 



- 1 



(10) 



which is valid for T^/t' > 1. Similar scaling arguments 
are used for the sum in Eq. lO which lead to 



le{x) 



sinTrfl 



fl 



„2e 



+ (1 - .x)2e + 2a;'' (1 - x) cos ttO 



(8) 

for < X < 1. The transition between the U shape 
behavior and the W shape behavior happens at Oc = 
0.5946.... The Lamperti PDF is shown in Figs. |2and|HI 
for the case r — Q, together with the numerical results. 



odd,ri<t' 



T^ 



it') 



/\l-9 



(11) 



an approximation which is vaHd for t' < T+. For t' > T+, 
N'^ — and X^Todd T<t' ^* ~ -^^- ^^ summary and after 
some rearrangements we obtain 



J 



CTA{t',r) 



^[0,T] j 1 ^ (l ^-^10,T]] 



T2 

[O.T] 



(1— r)I[u J,] 



1 - 



ire (l-r)I 



[0,T] 



r 



t' < ^^ 
t' > ^^ 



(12) 



Eq. lfT2|l yields the correlation function, however unlike 
standard ergodic theories the correlation function here is 
a random function since it depends on X^o^t]- The distri- 
bution of CrAit', T') is now easy to find using the chain 



rule, and Eqs. (jZIHlinil- In Figs. EElwe plot the PDF 
of CTA{t',T') (dashed curves) together with numerical 
simulations (diamonds) and find excellent agreement be- 
tween theory and simulation, for the cases where our ap- 




Figure 3: Same as Fig. [2| however now 6 = 0.8. If compared 
with the case 6 — 0.3, the distribution function exhibits a 
weaker non-ergodic behavior, namely for r = the distribu- 
tion function peaks on the ensemble average value of 1/2. 



proximations are expected to hold r < 1/2. We observe 
that unlike the r — case the PDF of the correlation 
function exhibit a non-symmetrical shape. To under- 
stand this note that trajectories with short but finite to- 
tal time in state on (r+ <C T) will have finite correlation 



functions when t' = 0. However when t' is increased the 
corresponding correlation functions will typically decay 
very fast to zero. On the other hand, correlation func- 
tions of trajectories with T+ ~ T don't change much 
when t' is increased (as long as t' ^ T+). This leads 
to the gradual nonuniform shift to the left, and "absorp- 
tion" into CTA{t',T') = 0, of the Lamperti distribution 
shape, and hence to non-symmetrical shape of the PDFs 
of correlation function, in general, whenever r ^ 0. 

Finally, we turn to the case T <^ t' . Since t' is large 
we use a decoupling approximation and write Eq. QJ as 



CTA{t',T') 



■ 1['d,T]1[t\T'\- 



(13) 



We distinguish between two types of trajectories, those 
in which no transition event occurs in the time interval 
[T, T'] and all other trajectories. Let Po[a, b) be the prob- 
ability of making no transition between time a and time 
6, also called the persistence probability [l^ . 



Po{a,h) 



sinTT^ 



a/b 



-^{i-xy'Ax 



(14) 



in the scaling limit. Using the Lamperti distribution for 
I[o,T] ) and probabilistic arguments with details left to the 
Appendix, we find the PDF of CTA{t' , T) 



J 



U.(t',T') {z) c. [1 - Po (T, T')] [1 - Po (f, T')] 



k {x) 



dx + M^ [;, (z) + 5 (z)] 1 -f Po (P, T') 
1 


[zleiz)^'^^] 
(15) 



Note that to derive Eq. ltT5|l we used the fact that I[o,t] 
and 2'[t',T'] are correlated. In Figs. I2I3I we plot these 
PDFs of CTA{t' ,T') (solid curves) together with numeri- 
cal simulations (diamonds) and find good agreement be- 
tween theory and simulation, for the cases where these 
approximations are expected to hold, r > 1/2. In the 
limit t'/T' ^ 1 Eq. ^ simplifies to 



Pc^,it',T'){z)--Mz)+6{z)]/2, 



(16) 



a result which is easily understood if one realizes that 
in this limit ^[t',T'] in Eq. (J13II is either or 1 with 
probabilities 1/2, and that the PDF of I[o,t] is Lamperti's 
PDF Eq. ©. 

To summarize, our work classifies the nonergodic prop- 
erties of photoemission intensity signal from NCs, and 
more generally Levy walks, and yields an analytical tool 



for the investigation of the nonergodic correlation func- 
tions. 
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Appendix A: DERIVATION OF EQ.(15) 

To calculate the PDF of CrAit', T') in Eq.(13) we use 
two steps: (i) calculate the PDF oiX\^t',T'] which statisti- 
cally depends on 1[o^t\ and then (ii) using the distribution 
of I[oT]) which is the Lamperti's PDF Eq.(8), calculate 
the PDF of CTA(i', Pa- 
using the persistence probability, we approximate the 
conditional PDF of Z^j/ j^/j for a given P[o,t] in the case 
T < i' by 



/i,,,,,, (^|I[o,T]) ^ [1 - ^0 (T, T')] Qi,,,,,, {z) + Po (T, T') [X[o,T]<^ (^ " 1) + (l ^ I[o,t])]^ {z)] 



(Al) 



where Qi ^, {z) is the PDF of I[t'^T'] conditioned that at least one transition occurs in \T,T']. In Eq. IjAljl we 



introduced the correlation between ^[t',T'] and Iio^t] through the dependence of the right hand side of the equation 
on 2^[o,T]- We assumed that in the case of no transitions in the time interval \T,T'], the probability of the interval 
[i',T'] to be all the time either on or off (the only possible choices) is linearly proportional to the value of 2r[o,T]- 
The persistence probability controls also the behavior of 



Qi 



[z) ^ [1 - Po {t\ T')] e (0 < z < 1) + Fo (t', T') 



5{z) + S{z-l) 



(A2) 



Briefly, we assumed that if a transition occurs in the 
interval [t',T'] the distribution oiX^t',T'] is uniform [i.e., 
9 (0 < z < 1) = 1 if the condition in the parenthesis is 
correct]. This is a crude approximation which is, how- 
ever, reasonable for our purposes (however when 9 ap- 
proaches 1, this approximation does not work). The delta 
functions in En. ljA2|l arise from two types of trajectories: 
If no transition occurs either 1[t' ,t'] = 1 (state on) or 



^[t',T': 



(state of f) with equal probability. 



Based on Eq.(13), the PDF of CxAit' , T') is 



PcTA{t',T'){z) ~ / te{x)fi^^,^,^ (J|x) 



dx 

X 



(A3) 



where we use the observation that ^[o.t] is distributed 
according to Lamperti distribution Eq.(8). Finally, from 
Eqs. I|A1IA2IA3|I we obtain Eq.(15). 
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